Weakly turbulent instability of anti-de Sitter space 
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We study the nonlinear evolution of a weakly perturbed anti-de Sitter (AdS) space by solving 
numerically the four-dimensional spherically symmetric Einstein-massless-scalar field equations with 
negative cosmological constant. Our results suggest that AdS space is unstable under arbitrarily 
small generic perturbations. We conjecture that this instability is triggered by a resonant mode 
mixing which gives rise to diffusion of energy from low to high frequencies. 
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Introduction. The past decade has witnessed growing in- 
terest in spacetimes which asymptotically approach anti- 
de Sitter space (AdS) which is the unique maximally sym- 
metric solution of the vacuum Einstein equations with 
negative cosmological constant. Despite this flurry of ac- 
tivity, motivated mainly by the AdS/CFT duality conjec- 
ture, the very basic question "Is AdS stable?" |H has been 
rarely raised (cf. (H for a notable exception) , let alone an- 
swered. This is in stark contrast to two other maximally 
symmetric solutions of the vacuum Einstein equations - 
Minkowski space (with zero cosmological constant) and 
de Sitter space (with positive cosmological constant) - 
which are known to be stable under small perturbations 
[3, The key feature of asymptotically AdS space- 
times, distinguishing them from asymptotically flat or 
dS spacetimes, is the presence of a timelike boundary at 
(spatial and null) infinity where suitable boundary con- 
ditions need to beprescribed in order to make the evolu- 
tion well-defined [5| . For no- flux boundary conditions one 
gets a Hamiltonian conservative system on an effectively 
bounded domain from which energy cannot escape 
(as for waves propagating inside a perfect cavity), hence 
the asymptotic stability of AdS space is precluded while 
the question of its stability, as we will see below, touches 
upon the KAM theory for partial differential equations. 

Model. In this Letter we report on numerical simulations 
which shed some new light on the problem of stability 
of AdS spacetime. As a simple model of asymptotically 
AdS dynamics we consider a self-gravitating spherically 
symmetric real massless scalar field in 3 -I- 1 dimensions 
whose evolution is governed by the Einstein-scalar system 
with negative cosmological constant A < 



Gap + Agais = SttG yda<j>di3<j> - ^gai3{d(pf 



For the metric we assume the ansatz 
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(1) 

(2) 



{-Ae-'^^dt^ + A-^dx^ + srn^x d^f) , (3) 



where P — —3/ A and dVP' is the standard metric on 
the round unit two-sphere. The ranges of dimensionless 



coordinates are — oo < t < oo (for the 'unwrapped' ver- 
sion of AdS considered in this paper) and < a; < 7r/2. 
We assume that A, 5, and </> are functions of {t^x). In- 
serting the ansatz ^ into the field equations (1-2) and 
introducing auxiliary variables <^ — cf)' and 11 = A~^e^(l) 
(hereafter overdots and primes denote derivatives with 
respect to t and x, respectively), we obtain a coupled 
quasilinear elliptic-hyperbolic system, consisting of the 
wave equation 



n 



tan^x 



{lecn^xAe-^'^)', (4) 



and the constraints (using units where 47rG = 1) 



A' = 



1 + 2 sitTx 



sm X cos X 



(1- A)-sinxcosxA($2+n2)^ (5) 



5' = - sin X cos x (<f>^ + 11^ 



(6) 



Note that the length scale £ drops out from the equations. 
The pure AdS solution is given by = 0, A = 1,(5 = 0. 
We want to solve the system (4-6) for smooth initial data 
with finite total mass. Smoothness at the origin implies 
that near x — Q we have the power series expansions 



<t>{t, x) - /o(<) + 0{x^), S{t, x) = 0{x^), 
A{t,x) = 1 + C'(a;2) , 



(7) 



where we used normalization 5(t, 0) = so that t is the 
proper time at the center. These expansions are uniquely 
determined by a free function fo{t). Smoothness at spa- 
tial infinity and finiteness of the total mass AI imply that 
near x — tt/2 we have (using p = tt /2 ~ x) 

^{t, X) = foo [t) p3 + O , 5{t, X) = Soo{t) + O , 

A{t, x) = l- 2Mp3 + o (p^) , (8) 

where M and free functions /oo(i), Soo{t) uniquely de- 
termine the power series expansions. It follows from 
that for smooth initial data there is no freedom in im- 
posing the boundary data (this is directly related to the 
fact that the corresponding linearized operator is essen- 
tially self-adjoint [Z]), hence the initial value problem is 
well-defined (see [8| for the rigorous proof) without the 
need of specifying boundary data at infinity. 
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Numerical results. We solved the system (4-6) numeri- 
cally using a fourth-order accurate finite-difference code. 
We used the method of lines and a 4th-order Runge- 
Kutta scheme to integrate the wave equation Q in time, 
where at each step the metric functions were updated 
by solving the hamiltonian constraint ^ and the slicing 
condition ^ . Preservation of the momentum constraint 
A + 2 sinx cos a: A^e~''$ n = was monitored to check 
the accuracy of the code. 

Solutions shown in Figs. 1 and 2 were generated from 
Gaussian-type initial data of the form 



$(0,a;) = 0, 



2e 

n(0, x) = — exp 

TT 



4 tan"' a; 



(9) 



with fixed width a = 1/16 and varying amplitude e. For 
such data the scalar field is well localized in space and 
propagates in time as a narrow wave packet. For large 
amplitudes the wave packet quickly collapses, which is 
signalled by the formation of an apparent horizon at a 
point XH where A{t, x) drops to zero. As the amplitude 
is decreased, the horizon radius xh decreases as well and 
goes to zero for some critical amplitude £o . This behavior 
is basically the same as in the asymptotically flat case, 
because for xh <C 7r/2 the influence of the AdS bound- 
ary is negligible. At criticality the A term becomes com- 
pletely irrelevant, hence the solution with amplitude 
asymptotes (locally, near the center) the discretely self- 
similar critical solution discovered by Choptuik in the 
corresponding model with A = [9[. For amplitudes 
slightly below eg the wave packet travels to inflnity, re- 
flects off the boundary, and collapses while approaching 
the center. Lowering gradually the amplitude we find 
the second critical value ei for which xh =0. As e keeps 
decreasing, this scenario repeats again and again, that is 
we obtain a decreasing sequence of critical amplitudes e„ 
for which the evolution, after making n reflections from 
the AdS boundary, locally asymptotes Choptuik's solu- 
tion. Specifically, we verified that in each small right 
neighborhood of En the horizon radius scales according 
to the power law xh{£) ^ (e — e„)^ with 7 ~ 0.37. Fig. 1 
shows that Xh{£) has the shape of the right continuous 
sawtooth curve with finite jumps at each e„. Notice that 
T(£„+i) — T(£„) « TT, where T(e) denotes the time of 
collapse. We stress that Xh is the radius of the first ap- 
parent horizon that forms on the t = const hypersurface; 
eventually all the matter falls into the black hole and 
the solution settles down to the Schwarzschild-AdS black 
hole with mass equal to the initial mass M (cf. [l^!)- It 
appears that lim„_^oo £n = 0, indicating that there is no 
threshold for black hole formation, however we did not 
determine precise values of £„ for n > 10 because the 
computational cost of bisection increases rapidly with n 
(since, in order to resolve the collapse, solutions have to 
be evolved for longer times on finer grids ). 

Let us mention that the analogous problem in 2 -I- 1 di- 
mensions was studied previously by Pretorius and Chop- 



X 




FIG. 1: Horizon radius vs amplitude for initial data The 
number of reflections off the AdS boundary before collapse 
varies from zero to nine (from right to left). 



tuik 11| who emphasized the challenges inherent in nu- 
merical simulations of AdS dynamics, however their anal- 
ysis was primarily focused on the threshold for black 
hole formation before any reflection off the AdS boundary 
takes place (as for our data with amplitude £0). 

In the following we consider the development of gen- 
eral (gaussian and other) small initial data, focusing at- 
tention on early and intermediate pre-collapse phases of 
evolution. We found that the Ricci scalar at the center, 
R{t,0) ^ ~2lP{t,0)/P-l2/e'^, can serve as a good indi- 
cator for the onset of instability. This quantity oscillates 
with frequency ~ 2 (as it takes time « tt for the wave 
packet to make the round trip from and back to the cen- 
ter) . An upper envelope of these oscillations is shown in 
Fig. 2a, where several clearly pronounced phases of evolu- 
tion can be distinguished. During the first phase the am- 
plitude remains approximately constant but after some 
time there begins a second phase of (roughly) exponen- 
tial growth, followed by subsequent phases of steeper and 
steeper growth, until finally the solution collapses. We 
find that the time of onset of the second phase scales as 
£^^ (sec Fig. 2b), which means that arbitrarily small per- 
turbations eventually start growing. Note that this be- 
havior is morally tantamount to instability of AdS space, 
regardless of what happens later, in particular whether 
the solution will collapse or not. In the remainder of this 
Letter we sketch a preliminary attempt to explain the 
mechanism of this instability in the framework of weakly 
nonlinear perturbation theory. 

Weakly nonlinear perturbations. We seek an approximate 
solution of the system (4-6) with initial data {(/), 0)|t=o = 
{ef{x),eg{x)) in the form 



2j + 



2i 



(10) 
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FIG. 2: (a) n^(t, 0) for solutions with initial data ((9| for four 
moderately small amplitudes. For clarity of the plot only 
the upper envelopes of rapid oscillations are depicted. After 
making between about fifty (for e = 6\/2) and five-hundred 
(for £ = 3) reflections, all solutions finally collapse, 
(b) The curves from the plot (a) after rescaling e^^n^(e^t, 0). 

where (j)i)t=o — (/(s^), 5(2;)) and all higher-order iter- 
ates have zero initial data. Inserting (jlOp into the system 
(4-6) and collecting terms of the same order in e, we ob- 
tain a hierarchy of linear equations which can be solved 
order-by-order. At the first order we obtain 

4>i+L(f>i^0, \-d^(tan^xd^) . (11) 

tan^x 

This equation is a particular case of the master equa- 
tion describing the evolution of linearized perturbations 
of AdS space, analyzed in great detail by Ishibashi and 
Wald [Tj. The Sturm-Liouville operator L is essentially 
self-adjoint on L^([0, •7r/2], tan^x da;). Below we denote 
the inner product on this Hilbert space by (/, 5) := 
j^^"^ f{x)g(x)ta.v?xdx. The eigenvalues and eigenfunc- 
tions of L are w| = (3 + 2jY (j = 0, 1, . . . ) and 
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ej{x) = dj cos^x 2-Fi(-J, 3 + j, -; sin^cc) , (12) 



1/2 

where dj = (16(j + + 1)It^)^ is the normalization 
factor ensuring that (ei,ej) — bij. The positivity of all 
the eigenvalues implies that AdS space is linearly stable. 
By elementary separation of variables, the solution of 
Eq. ()ll|) is given by the superposition of eigenmodes 

CXD 

(t, x) = flj cos(wjt -f /3j) Cj (x) , (13) 

where the amplitudes aj and phases jij are determined 
by the initial data. 

The back-reaction on the metric appears at the second 
order and can be readily integrated to yield 

COS'^X / • \ 

A2{t,x) = -. / Ui(t,2/)2-f(/,;(t,y)2 tan^ydy, 

sm X J Q \ / 

(14) 

62it,x) = - J (j)i{t,yf + (t>[{t,yf^ sinycosydy. 

(15) 

At the third order we get the inhomogeneous equation 

•jig +i03 = 5(01,^2,(52), (16) 

where S := -2(^2 + S2)4>i ~ (^2 + 62)^1 - (^2 + S'^Wi- 
Projecting Ea. ([T6)) on the basis ([T2|) we obtain an infi- 
nite set of decoupled forced harmonic oscillations for the 
generalized Fourier coefficients Cj{t) :— {(p^^ej) 

Cj + ^]cj = Sj (5, ej) . (17) 

Let I — {j G No : aj 7^ 0} be a set of indices of 
nonzero modes in the linearized solution ([T3| . A lengthy 
but straightforward calculation yields that each triad 
(ji: i2, J3) € such that ujj = ujj^ +iijj^—u)j^ gives rise to 
a resonant term in Sj (i.e., a term proportional to coswji 
or smujjt). Some of these resonances can be removed by 
a multiscale (or Poincare-Lindstedt) technique, but for 
general initial data there is no way to remove all the reso- 
nances and consequently we get secular terms which grow 
linearly in time and invalidate the perturbation expan- 
sion when e^t = 0(1). For example, for the single-mode 
initial data, {(f>,(/))t=o = £(eo(x),0), only 5*0 contains a 
resonant term and this term can be eliminated by intro- 
ducing a slow modulation of phase in (|13l) to yield the 
approximation 4>{t,x) ~ £C0s(3t + ^s^t) eo{x) which is 
uniformly valid up to at least t — 0{e~'^). In contrast, for 
the two-mode initial data, {4>,4>)t=Q = £(eo(x)-|-ei(x),0), 
the resonant terms in Sq and Si can be eliminated by 
suitable phase modulations but the resonant term in S2 
persists and produces a secular term in C2{t) ^ tsinlt. 

In accord with this analysis, we observe a striking dif- 
ference in the evolution of these two kinds of solutions. 
For (non-resonant) one-mode data the solution remains 
near the initial state during the entire simulation time, 
which indicates stability (although we cannot exclude a 
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metastable behavior with a very long hfetime). In con- 
trast, for (resonant) two-mode data we observe the on- 
set of exponential instability at time t = 0{e~^), as in 
the case of gaussian-type data Notice that for n- 
mode initial data the number of triads satisfying the 
diophantine resonance condition grows rapidly with n 
and hence, so does the number of unremovable secular 
terms. We believe that these secular terms are progen- 
itors of the higher-order resonant mode mixing which 
shifts the energy spectrum to higher frequencies. To 
quantify this effect, we introduce a spectral decomposi- 
tion of the total energy as follows. We define the projec- 
tions $j := (VA^.e'j) and Ilj := (a/AII, e^), and next, 
using the orthogonality relationships (e^-,e'j,) = ujjSjk, 
we express the total conserved energy as the Parseval 
sum M = i/p"/' {A^^ +AIl^)tan^xdx = ET=aEj{t), 



where Ei 



n 



'^j can be interpreted as the j- 



mode energy |13l | . Numerical evidence corroborating the 
conjectured mechanism of instability is shown in Fig. 3 
which depicts the transfer of energy to higher frequencies 
for a solution whose energy is initially deposited in the 
first two modes, i.e. Ej{0) « for j > 2. 

The weakly nonlinear multiscale perturbation analysis 
of the system (4-6) will be elaborated in detail elsewhere 
12l |. where we will also connect the nonlinear instability 
of generic quasiperiodic solutions of the linearized equa- 
tions to the violation of non-resonance conditions of the 
KAM theorem for PDEs and Arnold diffusion. 



tions. An important difference between wave equations 
defined on a fixed background and Einstein's equations 
should be stressed. For the former, the weakly turbulent 
behavior is compatible with global-in-time smooth evo- 
lution. In contrast, for Einstein's equations the transfer 
of energy to high frequencies cannot proceed indefinitely 
because concentration of energy on smaller and smaller 
scales inevitably leads to the formation of a black hole. 

Admittedly, the results presented above raise more 
questions than give answers. One of the key physical 
questions concerns the role of negative cosmological con- 
stant A in the observed phenomenon. Is it dynamical or 
kinematical? In other words, is an extra attractive force 
due to A essential in triggering gravitational collapse for 
arbitrarily small perturbations, or the only role of A is to 
confine the evolution in an effectively bounded domain? 

In this paper we studied perturbations of AdS^ for the 
Einstein-massless-scalar system in 3 -I- 1 dimensions, but 
we observed qualitatively the same behavior (in particu- 
lar, instability of AdS^) for the 4-f 1 dimensional vacuum 
Einstein equations within the cohomogeneity-two biaxial 
Bianchi IX ansatz of 16] . This result and its implications 
for the AdS/CFT conjecture will be discussed elsewhere. 
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FIG. 3: Evolution of the fraction of the total energy contained 
in the first (A: -I- 1) modes :— jj J2'j=o two-mode 

initial data <f){Q,x) — e ^^eo(a;) -I- ^ei(a;)j with e — 0.088. 

Conclusions. There is growing (numerical and theoreti- 
cal) evidence that the phenomenon of weak turbulence, 
i.e. the tendency of solutions to shift their energy from 
low to high frequencies, is common for (non-integrable) 
nonlinear wave equations on bounded domains, notably 
it has recently been proven for the nonlinear Schrodinger 
equation on torus 14,1^. To our knowledge, this paper 
is the first result in this direction for Einstein's equa- 
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